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ABSTRACT 

We study the evolution of Low Mass X-ray Binaries hosting a neutron star and of 
C^l , millisecond binary radio pulsars, with the help of numerical simulations that keep 

into account the detailed evolution of the companion star, of the binary system and 
of the neutron star. According to general relativity, when energy is released during 
accretion or due to magnetodipole radiation during the pulsar phase, the system loses 
gravitational mass. Moreover, the neutron star can collapse to a black hole if its mass 
exceeds a critical limit, that depends on the equation of state of ultradense matter 

■ and is typically ~ 2M Q . These facts have some interesting consequences: 1) In a 

millisecond radio pulsar the mass-energy is lost with a specific angular momentum 

I/"") , that is smaller than the specific angular momentum of the system, resulting in a 

■ positive contribution to the orbital period derivative. If this contribution is dominant 

and can be measured, we can extract information about the moment of inertia of 

. the neutron star, since the energy loss rate depends on it. Such a measurement can 

1 1 therefore help to put constraints on the equation of state of ultradense matter. 2)In 

Low Mass X-ray Binaries below the bifurcation period (~ 18 h), the neutron star 
survives the "period gap" only if its mass is smaller than the maximum non-rotating 
mass when the companion becomes fully convective and accretion pauses. Since in 
such evolutions ~ 0.8Mq can be accreted onto the neutron star, short period (P < 2h) 
millisecond X-ray pulsar like SAX J1808. 4-3658 can be formed only if either a large 
k> part of the accreting matter has been ejected from the system, or the equation of state 

of ultradense matter is very stiff. 3) In Low Mass X-ray binaries above the bifurcation 
period, the mass-energy loss lowers the mass transfer rate. As side effect, the inner core 
of the companion star becomes ~ 1% bigger than in a system with a non-collapsed 
primary. Due to this difference, the final orbital period of the system becomes 20% 
larger than what is obtained if the mass-energy loss effect is not taken into account. 

Key words: Stars: neutron - X-rays: binaries - binaries: close - pulsars:general - 
relativity - stars: individual: SAX J1808.4-3658 



1 INTRODUCTION 

Low-Mass X-Ray Binaries (LMXBs) are systems consisting 
of a neutron star (NS) with a relatively weak magnetic field 
(< 10 10 G) accreting from a low mass (~ IMq) companion 
star. When the companion star fills its Roche lobe, it trans- 
fers mass to the NS. The companion fills its Roche lobe 
either because it expands due to nuclear evolution or be- 
cause the lobe shrinks due to orbital angular momentum 
losses caused by gravitational radiation and magnetic brak- 
ing. The matter flowing from the inner lagrangian point to- 
wards the NS forms a Keplerian accretion disc around it. 
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The NS is spun up by the accreting matter, to an equilib- 
rium period that is roughly equal to the keplerian frequency 
at the inner rim of the accretion disc ( Ghosh 'fc Lamb 197911 . 
Once accretion ends, the NS can light up as a fastly ro- 
tating magnetodipole (radio pulsar): this is the so-called 
recycling scenario for the formation of millisecond pulsars 
( Bhattach arya fc van den Heuvel 1991| > . 

The secular evolution of LMXBs can follow two very 
different paths according to the evolutionary stage of the 
companion at the start of the mass transfer. If the orbital 
period of the binary at the beginning of the mass transfer 
is large, it begins when the companion evolves off the main 
sequence, the main driving mechanism for mass transfer is 
nuclear evolution and the system will evolve towards large 
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orbital periods. Such systems are said to be above the bi- 
furcation period (Tutukov et al. 1985, Ergma et al 1996). 
If the orbital period of the system at the onset of the mass 
transfer is small (e.g. it is below the bifurcation period), 
the companion is relatively unevolved, and the only impor- 
tant mechanism driving mass transfer are systemic angu- 
lar momentum losses (AML) due to magnetic braking and 
gravitational radiation. This type of evolution is similar to 
the classical evolution of cataclysmic binaries. The orbital 
period becomes shorter and shorter and the systems may 
experience a period gap when magnetic braking stops be- 
ing effective (when the secondary becomes fully convective) , 
and ultimately reaches a minimum period just before hy- 
drogen burning is extinguished (Paczynski & Sienkiewicz 
1981; Rappaport, Joss, & Webbink 1982). Beyond the pe- 
riod minimum (which depends on the evolutionary stage of 
the initial model), the donors radius begins increasing, fol- 
lowing the mass- radius relation for degenerate stars, and the 
orbital period will increase again, driven by gravitational ra- 
diation alone. Actually, the distribution of orbital periods of 
LMXBs, based on relatively few systems (Liu, van Paradijs 
& van den Heuvel, 2001), does not show a period gap as clear 
as that of CVs, but a more general limitation in the number 
of systems below ~ 4hr. It may well be that the distribution 
of initial periods of LMXBs favours the evolution of more 
evolved donors (Nelson & Rappaport 2003) and population 
synthesis results show that several LMXB systems may be 
descendants of intermediate mass secondary, after a phase of 
thermal timescale mass transfer (Podsiadlowski, Rappaport 
& Pfhal 2002), but the evolutionary path towards shorter 
periods is mainly followed by evolutions starting from un- 
evolved secondaries. For these systems, the donors have the 
same structure as the donors in CVs, and in principle there is 
no reason why they should not pass through a similar period 
gap, possibly having a smaller width, due to the fact that 
the accreting component is more massive (1.35 Mq) than 
the typical white dwarf in CVs (~ 0.6 Mq) (Podsiadlowski, 
Rappaport & Pfhal 2002). 

Studies of the evolution of LMXBs usually disregard the 
compact nature of the primary. This is not legitimate, since 
once the matter is transferred from the companion to the 
NS it releases mass-energy due to the strong gravitational 
binding of the NS, and therefore its gravitational mass (i.e. 
the "charge" of the gravitational force) decreases once it is 
accreted onto the NS. Thus even if the mass transfer is con- 
servative the total gravitational mass of a LMXB decreases 
when matter is transferred from the companion onto the 
NS: it has been suggested that this could have observable 
effects ( Alecian & Morsink 2004). Moreover, if enough mass 
is accreted onto the primary, it can collapse to a black hole 
(Lavagetto et al. 2004, hereafter Paper I). We should then 
keep into account both these effects if we want to study 
accurately the evolution of a LMXB. 

In this paper we will first of all introduce the evolu- 
tion equations for the binary system including the effects of 
general relativity (section 2). Then we will show how the 
mass-energy loss is potentially measurable in binary mil- 
lisecond radio pulsars (MSPs), and how this measurement 
can put constraints on the equation of state (EOS) of NSs 
(section 3). Moreover we will show in section 4 under which 
conditions the NS in systems below the bifurcation period 
can survive the "period gap". Finally in section 5 we will 



study the evolution of a system above the bifurcation pe- 
riod, showing how mass-energy loss alters the evolution of 
the system. 



2 THE EVOLUTION EQUATIONS 

We implemented a simulation code that includes the stellar 
evolution of the companion, the evolution of the binary sys- 
tem and the evolution of the NS under the effect of accretion 
at the same time in order to simulate accurately the evolu- 
tion of a LMXB. Our evolution code couples the routines 
of the ATON code (D'Antona, Mazzitelli and Ritter, 1988), 
updated with the physical inputs described in (Ventura et 
al. 1998), which accounts for the stellar evolution and the 
binary evolution of the system, and routines accounting for 
the evolution of the NS (which is considered to be fully rel- 
ativistic) under accretion . 

If the primary of the system is a NS, its gravitational 
mass will be given by the sum of the baryonic mass (i.e. 
the number of baryons N times the average bare mass of 
the baryons ms) and of the potential and kinetic energies 
divided by c 2 , that are non negligible since the gravitational 
binding energy is large for matter as dense as neutron star 
matter. We know that for any given equilibrium configura- 
tion of a NS we have on a general basis ( Bardccn 1970 ) 



M G = M G (M B ,J), 



(1) 



where we indicate the gravitational mass of the star with 
Mo, its baryonic mass with Mb and its intrinsic angular 
momentum with J. The accreted gravitational mass per unit 
time depends then both on the number of baryons accreted 
and on the accreted angular momentum IBardeen 
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where ujns is the NS spin frequency, c is the speed of light 
and $ is the energy needed to bring a unit mass from infinity 
to the pole of the star. 

Therefore equation becomes: 



M G = $M fl + ^ j. 



(4) 



When matter is transferred from the companion onto 
the NS, and the mass transfer is conservative, we have 



M B = -M c 



(5) 



where M c is the mass of the companion. It is useful to rewrite 
equation Q using equation as JAlecian fc Morsink "20041 



M G 



-(1-/3)M C 



(6) 



1 For a detailed description of the fully relativistic study of the 
response of the NS to the accretion of matter see Paper I. 
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where 



< p = 1 - $ 



UlNS J 



< 1. 



According to General Relativity the binding energy of the 
accreting matter results therefore in a mass defect. The grav- 
itational mass lost in accretion is released from the system 
(mainly as X-rays) and carries away a specific orbital angu- 
lar momentum that we assume to be equal to the specific 
orbital angular momentum of the NS: 



1 • 2rr / M c y 



(7) 



where P is the orbital period of the binary system and a is 
the orbital separation. Using the Kepler's law 
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where Mtot = M c + Mg , and the expression for the orbital 
angular momentum 

L = McM Q (-^-) 1/2 
we can rewrite equation J7J as 



(8) 
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where q — M c /Mg. The total variation of the orbital angualr 
momentum L will then be equal to the sum of the systemic 
orbital angular momentum losses L sys and of the angular 
momentum losses due to the relativistic mass defect Lg. We 
can therefore write: 
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where we made use of equation ©. If we now substitute 
equation © into equation I1UI . we can write the derivative 
of the orbital separation as: 



= 2 



2(1-9) 
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The last term on the right is due to the relativistic mass 
defect: it relevant only in compact systems where P is non- 
negligible. Using Kepler's Law, we can write for the evolu- 
tion of the orbital period: 
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(12) 



where the third term appears due to the relativistic mass 
deficit. This term is positive since — M c is positive: this 
means that as mass is transferred, an addictional positive 
contribution to the orbital period derivative is present in 
relativistic systems. One should keep in mind that the mass 
transfer rate M c depends upon the whole evolution of the 
binary system: in general it is a function of the nuclear evo- 
lution of the companion and of the orbital separation. In 



order to evaluate quantitatively the influence of the evolu- 
tion of the NS on the evolution of a system we must carry 
out detailed numerical simulations of the binary system due 
to the non-linearity of the equations. It is straightforward to 
extend these equations to the case of non-conservative mass 
transfer. 

In the next sections, we will show how these effects, 
together with the evolution of the compact object, can both 
have observable effects and alter the secular evolution of 
LMXBs. 



3 OBSERVABLE EFFECTS OF GENERAL 
RELATIVITY 

We may ask ourselves if the effects of relativity on the or- 
bital parameters of the binary system can be observable. 
In theory, any binary system about which we have enough 
information can reveal the effects we described in the pre- 
ceding section. Anyway, since we do not know much about 
most binary systems and, as various effects can overlap (see 
for example section 5), we have little chance to observe di- 
rectly these effects, which they only play a role in the secular 
evolution of the system. Alecian & Morsink (2004) argued 
that if the orbital period derivative of an accreting LMXB 
can be measured, it will allow to derive information on the 
structure of the NS, such as its mass and its gravitational 
binding energy. This effect can anyway be measurable only 
if the mass transfer is driven exclusively by the emission 
of gravitational radiation: there are too many uncertainties 
on the amount of angular momentum lost due to magnetic 
braking and on the nuclear evolution of the companion to 
allow to separate the relativistic effects in the orbital evolu- 
tion in such cases. This limits the range of systems that can 
be interesting to close LMXBs (with P<2 h) with a main 
sequence companion. In this situation the mass loss from 
the system will result in a potentially observable modifica- 
tion of the orbital period derivative. As Alecian & Morsink 
point out, however, uncertainties in the physics of binary 
evolution and of mass accretion make it difficult to separate 
this effect on the orbital period from others. Moreover, it 
is impossible to infer the mass accretion rate in a LMXB 
from its observed X-ray luminosity with the accuracy that 
is needed to extract information in their model. It is there- 
fore unlikely that this effect can be used to investigate the 
EOS of ultradense matter. 

On the other hand, when accretion onto the primary 
ends, the NS lights up as a radio pulsar, and it brakes down 
due to rotating magnetodipole emission. The NS loses some 
of its gravitational mass because it is radiating away energy. 
We can then write equation il l 21 in the form 
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where Pgw IP = SLqw/L is the orbital period derivative 
due to the emission of gravitational waves from the binary 
system, that we can write as tLorimer 2001t : 
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where e is the eccentricity of the system. It is then evident 
that when enegry is released from the MSP, the mass defect 
yelds a positive contribution to the orbital period derivative, 
opposite to the contribution of gravitational waves emission. 
Since the effect of gravitational mass loss is oc P, while the 
effect of gravitational waves emission is oc P _5//3 , the former 
effect will be dominant in systems with large enough orbital 
periods (say P > 6 h), while the latter will become more 
relevant in systems with short periods. 

The variation of the gravitational mass is equal to the 
spin-down energy of the pulsar, divided by c 2 . To a good 
approximation we can write: 



(15) 



where I is the moment of inertia of the NS. Esposito and 
Harrison (1975) readily noticed after the discovery of the 
first binary pulsar, PSR 1913+16, that the mass defect can 
alter the orbit of a binary pulsar. They found that its ef- 
fect on the orbital evolution of PSR 1913+16 was negligi- 
ble. Today we know many MSP, that have periods below 5 
ms and whose spin-down mass loss is larger (see equation 
I15|l . In these systems, the effect can be orders of magnitude 
stronger that the effect of gravitational waves. When this is 
the case, measuring the orbital period derivative of the bi- 
nary system can help to put strong constraints on the EOS 
of NSs. In many cases it is possible to measure both the 
spin frequency and its derivative in a radio pulsar with high 
precision. Moreover, the absence of mass transfer cuts away 
any uncertainty in the binary evolution model. The orbital 
period derivative depends then only on measured quantities 
(uims,^ns and e), on the masses of the two stars and on the 
moment of inertia of the neutron star (see equations 1131 1141 
and ll5H . We can get information on the two masses from the 
mass function 

f(M) = 



ML 



that is measurable in binary radio pulsars with very good 
precision. Using it, we can impose constraints on the mo- 
ment of inertia of the NS. Since the moment of inertia de- 
pends strongly on the EOS of the NS (ICook et al. 199411 . the 
detection of this effects will allow us to discriminate between 
various EOSs on a solid observational basis. 

Let us show this method with an example: suppose that 
we will observe a system with an orbital period P — 8 h, 
a spin period P s = 2 ms, P s = 3 x 10~ 19 , a mass function 
of 5 x 1O _3 M0 and that we have measured the orbital pe- 
riod derivative to be +2.5 x 10 -14 . In figure^we plot the 
values of the masses of the two stars that are compatible 
with this value of the orbital period derivative, in the hy- 
pothesis that the NS is governed by the pure neutron EOS 
by Pandharipande, named EOS A in the classic catalog by 
Arnett & Bowers (1977), or by the realistic hadronic EOS 
by Baldo, Bombaci and Burgio (1997), that we label as EOS 
BBB. The mass function imposes that the values of the two 
masses should be above the dotted line in figure □ The figure 
tells us that, given the mass function and the orbital period 
derivative we assumed, the pulsar cannot be described by 
EOS A. 

Now that we have shown how promising is this effect in 
principle, we may ask ourselves which is, between the known 



0.4 
0.35 



1 

— 


1 1 1 1 


1 

— 


— — ' 




BBB 










\a 




1 


1 1 1 1 


1 



1.3 



1.4 1.5 



1.6 1.7 

M G (M ) 



1.8 1.9 



Figure 1. Allowed values of the mass of the primary versus the 
mass of the companion (both in solar masses) for a system orbital 
period P = 8 h, a spin period P s = 2 ms, P s = 3 X 10 -19 and 
an orbital period derivative of 2.5 X 10 — 14 . The two lines in the 
figure are for NSs with EOS A (dashed line) and NSs with EOS 
BBB (solid line). The dotted line indicates the lower limit on the 
companion mass obtained if the mass function is 5 X 10~ 3 Mq. 



MSPs, the best candidate for detecting such an effect? The 
most promising object we found is PSR J0218+4232, which 
is constituted of a NS spinning at 2.3 ms in orbit with a white 
dwarf companion, with an orbital period of 2 days . The mass 
ratio is measured to be q = 0.13 ± 0.04 (Bassa et al. 20031. 

Due to its strong power output, the pulsar loses gravi- 
tational mass at the rate 4 X 10 _12 /45 M©/ yr, where J45 is 
the moment of inertia of the NS in units of 10 45 g cm 2 . The 
corresponding orbital period derivative, according to equa- 
tion 1131 . becomes P or b = 2.5 x 10~ 14 /45. The variation in 
the orbital period is derived from the measure of the peri- 
astron time delay. The effect of an orbital period derivative 
on the periastron arrival time is given by: 



AT, 



per — 0.5 — AT obs . 



(16) 



We find that we will need 117 /I^ 2 years of observation to 
detect a delay of 1 second, thus rendering impractical to 
measure such an effect in this pulsar, although period deriva- 
tives of the same order of magnitude have already been de- 
tected (see e.g. Nice et al. 2004), but in pulsars with short or- 
bital periods. Since the relativistic effect becomes dominant 
in pulsars with a large enough orbital period, we should seek 
for a pulsar with an higher spin-down energy - and therefore 
with an higher orbital period derivative for the same orbital 
period - in order to measure this effect in a shorter obser- 
vation time. Such a measurement is then quite unrealistic 
for presently known millisecond pulsars. However if we will 
find a pulsar with sufficiently high spin-down power it can 
become observable in a decade, becoming therefore feasible. 
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4 HOW ARE SHORT PERIOD LMXBS 
FORMED? 

Let us consider a binary system that is below the bifurca- 
tion period, and is constituted of a companion of 1M Q and 
of a primary with an initial gravitational mass of 1.35M0 
- a mass that appears to be the typical of isolated NSs 
jThorsett fc Ch akrabart y 1999) . The orbital period of the 
binary at the onset of the mass transfer is ~ 8 h. In the 
standard scenario, the system transfers mass because it loses 
angular momentum due to magnetic braking of the compan- 
ion, at a rate large enough to push the star well out of ther- 
mal equilibrium for periods shorter than ~ 4 h. When the 
companion becomes fully convective (typically at an orbital 
period ~ 3 h), magnetic braking is thought to stop, the com- 
panion recovers thermal equilibrium and the mass transfer 
ceases jSpruit fc Hitter 1983| l . The binary system will con- 
tinue to shrink due to gravitational waves emission, until it 
reaches an orbital period ~ 2 h, when mass transfer resumes. 
The evolution of this system in the period vs. accretion rate 
plane is shown in panel a of figure [5] We used the magnetic 
braking law by Verbunt & Zwaan (1981) with a braking in- 
dex of 0.5. In doing this first simulation, we neglected any 
evolution of the NS. 

What happens to the system if we keep into account the 
evolution of the primary? The fate of the NS in the period 
gap happens to be very interesting. An average ~ 0.8Mq 
have been accreted, and the NS is therefore rapidly spinning 
(see Paper I), and it will light up as a millisecond radio 
pulsar since accretion has stopped. Two general relativistic 
effects are important in this phase: 

(i) The loss of gravitational mass from the pulsar yields 
an additional positive contribution to the orbital period 
derivative (see the preceding section). This additional effect 
contrasts the shrinking of the orbit due to the gravitational 
waves emission, thus increasing the duration of the detached 
phase of the system. This increase can vary strongly depend- 
ing on the spin-down energy of the primary and on q (see 
equation 1131 . 

(ii) The silent collapse to a black hole if the pulsar 
is supramassive (i.e. its mass exceeds the maximum non- 
rotating mass), once it loses enough angular momentum. 

For the NS to survive the gap, the timescale of the collapse 
(T c ), must be larger than the timescale needed for the system 
to "cross" the period gap (T 9 ). Else, the NS will collapse to 
a black hole before the mass transfer resumes. This happens, 
for example, in our simulated system if the primary is a NS 
governed by EOS BBB, as can be seen in panel b of figure 

m 

The collapse time T c is defined by the equation 



Jin Jcrit 



Jdt 



(17) 



where Ji n is the angular momentum at the beginning of the 
detached phase, J cr it is the critical angular momentum be- 
low which the star collapses, and J is the angular momen- 
tum lost during the pulsar phase. Combining equation @ 
and the formula for the energy released by a magnetodipole 
rotator in general relativity (Rczzolla & Amhedov 2004), J 
is given by 
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Figure 2. Mass accretion rate as a function of the orbital period 
for a system below the bifurcation period. The companion is a 
population I star of IMq. We show the evolution of the system, 
not taking into account the evolution of the NS (dotted line) 
and the evolution of the system when we include the evolution 
of the primary: the companion transfers most of its mass due to 
magnetic braking, passing from P ~ 8 h to P ~ 3.5 h when 
the star becomes fully convective and mass transfer stops. The 
mass transfer then resumes once the system is close enough, at 
P ~ 2 h in the system with an unevolved primary. When we keep 
into acocunt the evolution of the primary, the system evolves 
similarly to the classical one, but when the companion becomes 
fully convective the NS is supramassive and therefore it collapses 
to a black hole before accretion can resume. As can be seen from 
the figure, no real difference is present in the evolution above 
the period gap, even if the relativistic mass defect is kept into 
account. For the sake of claryty, data were smoothed clean the 
numerical noise. 
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The gap time T g is defined by the equation 
Pdt 



(19) 



(20) 
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where AP gap is the amplitude of the period gap and P is 
defined in equation 11311 . In most situations AP gap ~ 1 h, 
and the mass of the companion is <0.25A/q. Integration 
of equations and JTJJ shows that T g > T c if the NS is 
supramassive and if the dipole magnetic field of the neutron 
star exceeds 10 26 G. This result holds for a vast range of 
EOSs from softer ones like EOS A to the stiffer ones like 
EOS L JArnett fc Bowers 1977L including recent, realistic 
EOSs like EOS FPS (|Lorenz, Ravehall fc Pethick 1993) , 
EOS APR lAkmal et al. 199811 and EOS BBB 
(Bald o~et al. 199"7l . as long as we assume that the 
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mass of the NS exceeds of at least O.O2M0 the maximum 
non-rotating mass. The result holds independently from 
the angular momentum of the NS at the onset of the 
pulsar phase. This means that LMXBs that host a NS and 
have a period <2 h, like all the millisecond X-ray pulsars 
known to date, cannot be supramassive if they have a 
non-negligible magnetic field. This result obviously holds 
only if the system is a NS-MS (main sequence) binary that 
has evolved through the period gap. 

For instance, the surface magnetic field of the first 
millisecond X-ray pulsar discovered, SAX J1808. 4-3658 
l |Wijnands fc van der Klis 1998^, has been estimated to be 
in the range (1-5) x 10 s G (IDi Salvo fc Burderi 20031 . If this 
system has a MS companion, so that it evolved from longer 
orbital periods, it cannot be supramassive, as it survived the 
period gap. The primary will not be a supramassive NS only 
if one of the following holds: 

(i) A relevant part of the accreting matter has been 
ejected from the system, and the mass transfer has therefore 
been non-conservative for most of the binary evolution. 

(ii) The maximum non-rotating mass of the NS is very 
high, > 2M Q (i.e. the EOS of the NS is very stiff). 



5 SECULAR EVOLUTION OF SYSTEMS 
ABOVE THE BIFURCATION PERIOD 

In the preceding section, we have shown how the evolution of 
the compact object can alter the evolution of a system below 
the bifurcation period. Now we will show how the evolution 
of the NS can alter the secular evolution of a system above 
the period gap. In such a system, the mass transfer is driven 
by the nuclear evolution of the companion. The evolution of 
such wide systems is well understood (We bbink et, al. 1983). 
and is thought to explain well the formation of binary MSPs 
with large orbital periods 

In order to show how this canonical evolution is altered 
when we keep into account the evolution of the primary, 
we simulated a binary system whose companion is a 1.1 Mq 
population I star. The initial mass of the primary is chosen 
again to be 1.35 Mq. When the mass transfer begins the 
orbital period of the system is 11 d. 

First of all, we carried out a simulation in which the 
we disregarded the evolution of the primary, for which we 
assumed Mg = Mb- This system is indicated in the fol- 
lowing as system A. We simulated also a system (named B 
in the following) with the same companion star, but this 
time taking into account the evolution of the primary. The 
NS was assumed to have a low surface magnetic field of 
10 8 G, and the EOS has been fixed to be EOS BBB we al- 
ready introduced in section 3. In paper I we showed that 
a weakly magnetized NS is easily spun up to periods well 
below one millisecond, uncomfortably lower than the mini- 
mum observed period for a NS, 1.56 ms (Backer et al. 1982 ). 
There is mounting evidence that some mechanism, whose 
nature is still unclear, prevents NSs to spin faster than 700 
Hz jChakrabar ty et al. 2003] l . We simulated therefore a sys- 
tem identical to system B, but we kept the spin frequency 
artificially below 700 Hz. We will refer to it as system C. 

As we said before, the evolution of the mass of a NS 
differs from that of a non-collapsed star because of the dis- 
crepancy between the gravitational and the baryonic masses 



in NSs. According to equation Q the gravitational mass of 
the NS will be smaller if the spin frequency and the angular 
momentum of the star are smaller for a given baryonic mass. 
Therefore, we expect that the mass of the NS in system C 
will be smaller than the one of the NS in system B, and both 
will smaller than that of system A. From our simulations we 
find in fact that the gravitational mass of the NS at the end 
of the evolution is 1.96 Mq in system A, 1.88 M (4% less) 
in system B, and 1.79 Mq (9% less) in system C. The higher 
mass defect in system C is due both to the smaller angular 
momentum of the star (that is almost irrelevant, since in 
a neutron star the rotational energy is usually ~ 0.1 times 
the binding energy) and to the fact that when the NS spins 
slower it is more compact and therefore its binding energy 
is also higher. 

It is interesting to note that the evolution of the NS 
has other effects on the evolution of the binary system. We 
know that the final state of such systems depends substan- 
tially only on the mass of the inner core of the companion 
star IWebbink et al. 198311 . Looking at equation (1121 we see 
that the relativistic effect is dominant during the first phases 
of accretion, when q ~ 1: when matter is transferred, a rel- 
ativistic system enlarges faster than a non-relativistic one. 
Consequently, the Roche Lobe of the companion is larger, 
and the mass accretion rate will become smaller. As a conse- 
quence, the relativistic systems have a longer time to evolve, 
and the mass of the core becomes slightly bigger: at the 
end of accretion, the companion in system A has a core 
of 0.305Mq. In system B the companion has a core that 
is 0.003Mq heavier than in system A, while in system C, 
the core is 0.006Mq heavier. As a consequence (see again 
Webbink et al. 1983), the orbital period at the end of the 
accretion process is 105 d in system A, 114.5 d (~ 9% larger) 
in system B and 124 d (~ 18% larger) in system C (as can 
be seen in figure [3J. 

This means that the net effect on the orbital period 
evolution of the system is stronger than the effect on the 
mass alone, but this is due to the small changes in the final 
core mass and not directly to the third term on the right 
of equation 111 2> . Small changes in the mass of the core also 
account for significant variations in the orbital period at the 
end of mass transfer: this means that, although the effects 
we describe can change the evolution of a system from given 
initial conditions, it does not alter the scenario of the evo- 
lution of such systems. 



6 CONCLUSIONS 

In this paper we showed that the presence of a NS, that can 
be described using general relativity, can have a big impact 
on the evolution of the binary system. 

First of all we searched for potentially observable ef- 
fects of the relativistic nature of the primary. We noticed 
that when the NS releases energy without accreting during 
the pulsar phase, it will lose gravitational mass. Therefore, 
a positive contribution to the orbital period derivative can 
dominate in certain situations over the negative contribution 
due to gravitational waves emission, resulting in an overall 
positive orbital period derivative. We showed that the mea- 
surement of such a period derivative in a binary millisecond 
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Figure 3. Orbital period (in days) of the binary system as a 
function of the mass of the companion (in solar masses). Letters 
indicate the evolutionary tracks for systems A,B and C. The evo- 
lution is followed only until the companion collapses to a white 
dwarf. The x axis is inverted, so that the evolution of the systems 
goes from left to right. 



pulsar can allow to put constraints on the EOS of ultradense 
matter on a solid observational basis. 

Then we concentrated on the effects that the secular 
evolution of the primary can have on the evolution of bi- 
nary systems, both below and above the bifurcation period. 
In systems starting below the bifurcation period a necessary 
condition for a NS to survive the period gap without collaps- 
ing to a black hole is to be non-supramassive. This means 
that the NS in LMXBs with short periods, like the millisec- 
ond X-ray pulsar SAX J1808. 4-3658, cannot be supramas- 
sive if they evolved from larger periods. This implies that 
either the mass transfer is highly non conservative or the 
EOS of ultradense matter is stiff. 

In systems starting above the bifurcation period, we 
showed how the relatively small effect of General Relativ- 
ity on the orbital evolution can alter the evolution of the 
companion; this means that the total effect can be non- 
negligible. However these effects do not change the proposed 
scenario for the evolution of large-period systems in any sig- 
nificant way. 

In all, the relativistic mass defect has some effect on the 
evolution of binary systems, but it is negligible if compared 
with the uncertainties that we still have on the other effects 
present in the standard theory of the binary evolution. For 
example, various laws have been proposed for the magnetic 
braking mechanism, and the intensity of the effect has been 
questioned (see e.g. Ivanova & Taam, 2003). However, the 
evolution of the primary - and the resulting mass defect - 
can be important in studying the evolution of the binary 
system, both because it can give rise to an observable effect 
in MSPs and because it puts constraints on the evolution of 
systems below the bifurcation period, and therefore cannot 
be disregarded when doing evolutionary studies. 
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